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Q^. Abstract We present a formulation of the Hamiltonian variational method for QED which 

<^ ! enables the derivation of relativistic few-fermion wave equation that can account, at least 

Q \ in principle, for interactions to any order of the coupling constant. We derive a relativistic 

T-H ■ two-fermion wave equation using this approach. The interaction kernel of the equation is 

' shown to be the generalized invariant M. matrix including all orders of Feynman diagrams. 

^ ■ The result is obtained rigorously from the underlying QFT for arbitrary mass ratio of the 

^ ! two fermions. Our approach is based on three key points: a reformulation of QED, the 

Q \ variational method, and adiabatic hypothesis. As an application we calculate the one-loop 

^ I contribution of radiative corrections to the two-fermion binding energy for singlet states with 

\^ • arbitrary principal quantum number n, and i = J = . Our calculations are carried out in 

O ■ the explicitly covariant Feynman gauge. 

^: 
Oh: 

^.' 1. Introduction 

_>: 

k> ! The description of two-particle states (particularly bound states) in QFT such as quan- 

H I tum electrodynamics, including relativistic effects, is an important problem. It is well known 

■ - - ' that in the nonrelativistic limit this problem has been solved on the basis of the Shrodinger 

equation. In the relativistic case we have the Dirac equation, which describes a one-fermion 
particle system only. However the implementation of this equation for two-particle or multi- 
particle system meets significant difficulties. It requires a relativistic manybody approach, 
which is provided by QFT. However, the QFT implementation to relativistic few-particle 
bound states is difficult. The usual method for treating two-body bound states in QFT is 
by means of the Bethe-Salpeter equation. It has a number of shortcomings, among them, 
a perturbative treatment of interactions, which is unsuitable for strongly coupled system. 
Other difficulties are the appearance of relative time coordinates, which make it problematic 
to interpret the Bethe-Salpeter amplitude as a traditional wave function. There are various 
ways of avoiding the problem, such as, e.g., single-time reductions, but these are generally 
non-rigorous approximations to the original problem in QFT. 

Recently we applied the variational Hamiltonian formalism to the problem of two fermions 



interacting through the electromagnetic field [1,2]. We derived two-fermion relativistic wave 
equations, classified all bound states and calculated the relativistic energy spectrum to fourth 
order in the fine structure constant. The results were obtained for arbitrary mass ratio of 
the particles. As a consequence, the fourth-order hyperfine splitting formulas were derived 
for arbitrary quantum numbers and mass ratio. These formulas are in agreement with 
experimental results for hydrogen and muonium, as well as with previous calculations. 

In the present paper we develop the variational Hamiltonian formalism for two-fermion 
systems to allow for the inclusion of effects of higher order, such as vertex corrections, vac- 
uum polarization, and two-photon exchange. This can be achieved by using the variational 
principle in combination with the adiabatic hypothesis. The two-fermion wave equations 
are derived rigorously from underlying QFT for arbitrary mass ratio. The solution of these 
equations allows one to calculate energy corrections to all orders on the basis of a two-body 
approach. The method presented here can be easily generalized to three- and multi-particle 
systems. It is important to note that the variational Hamiltonian formalism is applicable to 
strongly coupled systems for which perturbation theory may be unreliable. 

2. Reformulated QED 

It has been pointed out in previous publications that various models in Quantum Field 
Theory (QFT), including QED, can be reformulated, with the help of mediating- field Green's 
functions, into a form that is particularly convenient for the investigation of bound-state 
problems and variational calculations [3,4]. This approach was applied to the study of 
relativistic two-body eigenstates in the scalar Yukawa (Wick-Cutkosky) theory [5,6,7]. We 
shall implement such an approach to two-fermion states in QED in this paper. 

The Lagrangian for two fermionic fields interacting electromagnetically is {h = c = 1) 

C = ^{x) [tYdf, - mi - Q,YA^{x)) ^{x) + 0(a;) (^7^9^ - m^ - Q2l^A^{x)) <j){x) 

-i (d^A^ix) - dpA^ix)) {d^A^ix) - d^A^ix)) . (1) 

The corresponding Euler-Lagrange equations of motion are the coupled Dirac-Maxwell equa- 
tions, 

(,^/^5^ _ m,) ^{x) = QirA^{x)^P{x), (2) 

{tYd^ - m2) <p{x) = g27^^M(^)0(^)> (3) 

and 

9^9M-(x) - d^d,A^^{x) = fix), (4) 

where 

fix) = Q,'^ix)f^ix) + Q2^ix)f(j)ix). (5) 

Equations (2)-(4) can be decoupled in part by using the well-known formal solution [6-10] 
of the Maxwell equation (4), namely 

A^ix) = Alix) + [ d'x'D,,ix - x')fix'), (6) 



where D ^^{x — x') is a Green's function (or photon propagator in QFT terminology), defined 

by 

9,9"D^,(x - x') - d^d'^D^^ix - x') = g^j\x - x'), (7) 

and A^Jx) is a solution of the homogeneous (or "free field") equation (4) with j''(x) = 0. 

Equation (7) does not define the covariant Green's function D^^{x — x') uniquely. One 
can always add a solution of the homogeneous equation (eq. (7) with g^,^ -^ 0). This allows 
for a certain freedom in the choice of D^^,, as is discussed in standard texts (e.g. ref. [8,9]). 
In practice, the solution of eq. (7), like that of eq. (4), requires a choice of gauge. However, 
we do not need to specify one at this stage. 

Substitution of the formal solution (6) into equations (2) and (3) yields the "partly 
reduced" equations, 

{tYd^-m,)ij{x) = Q^^^(aI{x)+ f d'x'D^,{x-x')f{x')]iP{x), (8) 

(ij^d,-m2)(l>{x) = Q2J^(aI{x)+ fd'x'D,,{x-x')f{x')\(l){x). (9) 

These are coupled nonlinear Dirac equations. To our knowledge no exact (analytic or nu- 
meric) solution of equations (8)- (9) for classical fields have been reported in the literature, 
though approximate (perturbative) solutions have been discussed by various authors, partic- 
ularly Barut and his co-workers (see ref. [10,11] and citations therein). However, our interest 
here is in the quantized field theory. 

The partially reduced equations (8)-(9) are derivable from the stationary action principle 



SS [ip, (t)]=5 j (TxCr = (10) 

with the Lagrangian density 

jCr = ^{x) {t^d^ - mi - QiYAlix)) ^{x) + 0(x) {t^d^, - m^ - Q2l''Al{x)) <j){x) 

-\ld'x'f{x')D,,{x - x')f{x) (11) 

provided that the Green's function is symmetric in the sense that 

D^^{x -x') = D^^{x' -x), D^^{x -x') = D^^{x -x'). (12) 

The interaction part of (11) has a somewhat modified structure from that of the usual 
formulation of QED. Thus, there are two interaction terms. The last term of (11) is a 
"current-current" interaction which contains the photon propagator sandwiched between 
the fermionic currents. We shall use this modified formulation together with a variational 
approach to obtain relativistic two-fermion equations, and will study their bound-state so- 
lutions. 



We consider the quantized theory in the Hamiltonian equal-time formahsm. To this end 
we write down the Hamiltonian density corresponding to the Lagrangian (11), namely 



Hr = Ho+Hi, Hi = Tii^ + TY/j, (13) 



where 



Hq = ilj\x){-i'a -V + mif3)ilj{x) + (f)'^{x){~i'a -V + m2f3)(l){x), (14) 

Hi, = ^fd'x'r{x')D^,{x-x')f{x), (15) 

Hi, = QMx)rA'^{x)i,{x) + Q2^{x)YAl{x)<P{x), (16) 

and where we have suppressed the kinetic-energy term of the free photon field. We quantize 
the theory by imposing equal-time anticommutation rules for the fermion fields, namely 

^Jx,t),V^|3(y,t)} = {0jx,t),0l3(y,t)} =5,^53(x-y), (17) 

and all others vanish. In addition, there are the usual commutation rules for the Aq field, 
and the commutation of the Aq field operators with the ip and field operators. 

To specify our notation, we quote the Fourier decomposition of the field operators, namely 



with Pi = Pi = [ujp, p), Up = \/ml + p2 and 

'^(^) = E / 7^ (^) [^P^^ (P' ^) ^"^'"^ + ^IsV (P, s) e^^--] , (19) 



with P2 = P2 = (^P5 P)) ^p = V ^^2 + P^- Note that the mass-?T2i free-particle Dirac spinors 
u (p, s), V (p, s), which satisfy (7^Pi^ — rrii) u (p, s) = and (7^Pi^ + mi) v (p, s) = 0, are 
normalized such that 

m'^ (p, s) M (p, a) = v'' {p,s)v{p,a) = ^-Ssa, (20) 

mi 

u'' {p, s) V {p, a) = t;'f (p, s) M (p, 0-) = 0. 

Analogous properties apply to the mass-'m2 spinors U, V. The creation and annihilation 
operators 6^, 6 of the (free) particles of mass mi, and S, d for the corresponding antiparticles, 
satisfy the usual anticommutation relations. The non-vanishing ones are 

{6p„ h\^] = {rfp„ 44 = 5,J' (p - q) . (21) 

Again, the analogous properties apply to the mass-'m2 operators B\ B, D\ D. As a concrete 
example, we can think of the mass-mi particles as electrons, and the mass-'m2 particles as 
muons, though any pairs of charged fermions could be considered. 



3. Variational principle and adiabatic hypothesis 

The Hamiltonian formalism of QFT is based on the covariant eigenvalue equation 

P^ 1^) = P^ \i,) , (22) 

where P^ = iH,Pj is the energy-momentum operator, P^ = (E, P) is the corresponding 

energy-momentum eigenvalue, with E'^ — P^=M^, where M is the invariant mass of the sys- 
tem. The /i = component of (22) cannot be solved for realistic theories, and approximation 
methods must be used. We shall develop an approximation method in this section, based 
on the variational principle 

6[{^P,\H,it)-E\^p,),^,^)=0, (23) 

which is equivalent to the fi = component of (22), so to this end, we shall consider the 
expectation value. The subscript h is used to stress that the equation is written in the 
Heisenberg representation, which relates to the Schrodinger picture as \^f^) = e*^^(*~*") li's)- 
It is well known that variational solutions are only as good as the trial states that are 
used. Thus, it is important that the trial states possess as many features of the exact solution 
as possible. For our purpose it is convenient to rewrite equation (23) in the interaction 
representation 

6[{^,{t)\H.{t)-E\^,{t)),^,^^=0, (24) 

where the operator Hi (t) = /Jqi (t) + Hji (t) consists on two parts, which represent the 
free-field and interacting-field Hamiltonians respectively. Note that Hoi (t) = Hoi is time- 
independent. 

Using the time-evolution operator U we can express states in the form 

l^PM = U{tM) l^.(ti)) , mt)\ = {i^,{t,)\ U\tM), (25) 

The time-evolution operator U satisfies the familiar differential equation, 

Hn{t)U{t,to)=t-Uit,to), (26) 

with the boundary conditions U{to,to) = I, where / is the unit operator. Equation (26) has 
the iterative solution 

U (t, to) = f/(^) (t, to) + f/('^ (t, to) + f/(^) (t, to) + ..., (27) 

where 

^^'^(t,to) = /, (28) 

t/(2) (t, to) = -t [ dti f d^x.Hn (xi) , (29) 



f/(3) (t, to) = -^^ / dti j dh j d'^id^^2T (Un (xi) Hn {X2)) , (30) 

and Hj is defined by (13), (15) and (16). 
Substitution of (25) into (24) yields 

6(^{Ut2)\UHt,t2) (i/.(t)-i5)f/(t,ti)|^,(ti)),^,J =0. (31) 

Using tlie properties of tlie U operator, W(t,t2) = U{t2,t), and Ui{t2,i)Uj{t,ti) = Uj{t2,ti) 
we obtain for t2> ti 

5[{i^,{t2)\T{u{t2M) [H,,-Ey) +T{u{t2M)HH{t)) |^,(ti))^^^^^) = 0, (32) 

where T is the time-ordering operator. The first term in (32) is 

(^,(t2)|T(t/(t2,ti) [H,^-E)) \i^Ati))t=u, = {i^Ah)\Hu-E\i^,{h)), (33) 

since ifoi is a time-independent operator, and {'ipi(t2)\ U(t2,ti) = {ip^(ti)\. The second term 
in (32) can be rearranged as follows, 

{Ut2)\T{uit2,h)Hnit)) |^,(ti)),^,, = {^p,ih)\U-\t2,h)T {uit2,h)Hnit)) |^,(ti)),=,„ 

(34) 
If \ipi(t2)) and \ipi(ti)) are stationary states, they must coincide to within an arbitrary phase 
factor e*^, which does not affect any physical values. Hence 

f/(t2, ti) Mh)) = |^,(t2)) = e'' Mti)) , (35) 

where 

e'' = {^P,ih)\Uit2,h)Mh)). (36) 

Thus 

{Uti)\U{t2,h)Mti)) 
and we obtain 

/^ ^ X {Uti)\T(uit2,h)Hnit))Mt,)),^,^^ 

mt2)\T{u{t2,h)H,^{t)) mh)\^,, = nuMuu ,M u^^ ' ^^^^ 

It is known [12], [13] that the numerator of (38) can be written as a product of two factors, 
which correspond to all connected and disconnected diagrams. The factor, which represents 
the disconnected diagrams, is identical with the denominator {ilj^{ti)\ U{t2,ti) \ilj^{ti)), and 
therefore equation (32) takes the final form, 

5 ({Utl)\Ho^- EMt,)) + {Uti)\T (uit2,t,)Hnit)) |^,(ti)),= J = 0, (39) 



which involves connected diagrams only. This equation is just another form of (24) and does 
not contain anything new with respect to it. However, equation (39) is convenient for the 
implementation of adiabatic switching [12], [13], in which case we put ti = — oo. We assume 
that the Hamiltonian Hn (t) depends on an interaction constant a, hence we can switch this 
interaction on by the following anzatz 

a{t) = a ((1 - v) e-"l*-*»l + r]) , (40) 

where e and rj are small positive quantities. At time t = — cxo the coupling constant a (— oo) = 
at] is very small ("bait" interaction), and becomes the full interaction at t = to. This means 
that at t = — cxD the state |'?/'j(— oo)) is an eigenstate of a weakly interacting system, and by 
process of adiabatic switching, it becomes a stationary state \ipi(to)) = U(to, — oo) \tpi{—oo)) 
of the fully interacting system. This is the content of the adiabatic hypothesis. It is evident 
that the state |-?/'j(— oo)), in a sense, contains only minimal information about the interaction 
in the system, however it should include the main features of the "bait" interaction such as 
kinematical properties and symmetry of the system. The adiabatic hypothesis permits the 
use of simple trial states that take into account the key physical properties of the system. 

4. Trial state and two-fermion system 

The simplest Fock-space trial state, that can be written down in the rest frame (P = 0, 
E = M), for a two-fermion system is 

S1S2 

where b'^ and D_pg^ are creation operators for fermions of masses mi and m2 respectively, 
and |0) is the trial vacuum state such that bpsi |0) = -Dpsj |0) = 0. The -F^isa ^^^ four 
adjustable functions, that are chosen so that (41) is, as we show in section 5, an eigenstate 
of the relativistic total angular momentum operator, its projection, and parity. We shall use 
this trial state to implement equation (39), which we write in the form (ti = — oo, t2 = +oo, 
and to = 0) 

5 [{tPt„jHo, - E \^p,„J + {^Pt„ai\T (f/(+oo, -oo)Hn (t)) \^Pt„ai)t=o) = 0- (42) 

In practice U{+oo, — oo) is represented by the series (27), however any finite number of terms 
of the series may be included. In lowest order {U = I) equation (42) reduces to 

6 ((^,„,,| Ho,-E + Hu (t) \i'tr^al)t=o) = 0- (43) 

This equation was used in earlier works [1], [2] (with the trial state (41) and the interacting 
Hamiltonian Hu = Hj^i (15)) to obtain bound-state energy spectra to order a^ for all states 
of positronium- and muonium-like systems. In this work we shall go beyond this order. 



Note that only the term i7/^, eq. (15), of the interaction part of the Hamihonian con- 
tributes in eq. (43), since {iptriail -^^2* (^) li^triai) = with the choice (41) of trial state. That 
is, the simple trial state (41) does not sample the Hi^i interaction term, which means that 
processes that involve the emission or absorbtion of radiation (i.e., physical photons) are 
not accomodated. To do so would require a modification of the trial atate (41). However, 
this will not be considered in the present paper, which thereby deals with equations for pure 
bound states (or elastic scattering). 

The matrix elements (42) to implement the variational principle needed are 

{^P,„J :H,-E: \^P,„J =J2 J d'pF:^Jp)F,,s,{p) (u, + %-M) (44) 



SlS2 



and 



{,lj,„J : T (f/(+oo, -^)Hu (t)) : \^^r^al)t=o (45) 

(^,„,,| :y"d^xT((t/W(+oo,-oo)+t/(^)(+oo,-oo) + ...)7^,,(x)(5(t)) : |^,„J 



qiq2mira2 sr-^ f d pd q 



Yl / / ^ nn ^^*i^2 (P)^^i'^2 (q) i-i) Ms,s2a^a2 (P, q) , 



aia2SiS2~ V 'P^q^^P^^q 



(27r) 
where A^siS2o-io-2 (P) q) is the generalized invariant Ai -matrix 

Ms,s2.,a2 (p, q) = -^£^1.2 (p> q) + -^iiUia2 (p> q) + ••' (46) 

which includes effects in all orders of the interaction, and where the sum contains all Feynman 
diagrams, reducible and irreducible. Some details of the calculations of (45) to order Ai^'^^ 
are presented in Appendix A. 

The variational principle (42) leads to the following equation 

^ J £p (iu, + n,-E) F,,,,{p)6F:^Jp) (47) 



SlS2 



iTLim2 Y^ f d^pd^fi 



(2 



n 



\3 



I n n -^'^^"^(q) ("^) •^^i^2-i<x2 (P, q) ^Krs2^V) = 0. 

„ \/^n^o"n"o 



(Jlcr2SiS2 



We now discuss the structure of the Ai- matrix on the basis of the expansion (27) for 
U, and restrict our consideration to one-loop level only. In lowest order, U = U^^^ = I, we 
obtain 

-Mi'.U.2(P>q) = A^:f.V,., (P, q) _ (4J 

= -u (p, si) i-iQiY) u (q, CTi) iDf,^{p - q)V (-q, (73) {-iQ2Y) V (-p, S2) 

where A^s^fjo-io-a (P' q) is the usual invariant matrix element corresponding to the one-photon 
exchange Feynman diagram obtained and considered in [2]. This expression involves the 



Fourier transform of the Green's function (eq. (7)), namely 

D,Ax -^') = l -0y.D,Me~^'^--r (49) 

The Green's function D^y{p — q) consists of two parts 

D^^u{p ^ ^) = 2 (^M^(^i ~~ ^i) + ^^^M2 - P2)) , (50) 

where D^ij{pi — qi) and -D^jy(g2 — P2) are the Green's functions of particles with masses mi 
and 1712 respectively. Note that the four- vectors pi, p2, qi and g2 are defined as in formulae 
(18), (19). For a fermion-antifermion system like positronium we obtain [1] the additional 
virtual-annihilation term (Qi = Q2 = ^) for ■M.siS2aia2- 

A<s™2<xia2 (P, q) = M (p, si) i-ie-f^) V (-P, S2) iD^^ (up) v (-q, ^2) {-ie^) u (q, CTi) . (51) 

The terms of next order in the coupling require the inclusion of the operator f/*^^\ eq. (29), 
and can be written in the form 

K A ( 2 ) ^W ^^ / K A vac ' \ \ A VGV ' \ K A "^^^^^ a * j TfldSS - \ K A 2'Y / r' c\\ 

-^SlS20-l(T2 = 2-^ I -^31520-10-2 + -^SlS20102 +-^SlS2010-2 + -MsiS2O10-2 ) + -^SiS20lO-25 V^^) 

i=l,2 

where the index i = 1,2 corresponds to masses mi and m,2 respectively. One- loop level ra- 
diative corrections include second-order vacuum polarization Ai'""-'^, vertex corrections Ai'"'^^, 
and the fermion self-energy y\4™''*^. Their explicit form is given by 

-Miilroio2(p>q) (53) 

= -u (p. Si) (-iQi7^) u (q, ai) iD^^ {p - q) W"" {pi - qi) X 
xiD^p {p-q)V (-q, ^2) {-iQ2l^) V (-P, S2) , 

A^iilro7o2 (p> q) _ (54) 

= -u (p, Si) {-iQi) A" {pi, qi) u (q, 0-1) iDa/3 ip-q)V (-q, 0-2) {-iQ2l'^) V (-p, S2) , 

A^iifcTol (p, q) _ (55) 

= -u (p, Si) S (pi) S^ (pi) (-iQi7") u (q, ai) iD^p {p - q)V (-q, (72) {-iQ2l'^) V (-p, S2) , 

A^g&T^I (p, q) _ (56) 

= -u (p, Si) (-i(5i7") -^v, (pi) S (pi) u (q, (Ti) iD^^ (p - g) V (-q, (J2) {-iQ2l^) V (-p, S2) , 

where the standard definitions of the 11, A and S functions apply (we display functions and 
operators for the ip field only), namely 

n'^" (Pi - gi) = -tQl [ -^Tr [YS^ {k + pi- qi) rS^ (k)] , (57) 

J (27r) 

9 



for the vacuum polarization function, 

d*k 



A" (pi, gi) = -iQl I j-i^YS^ {k + p,) rS^ {k + gi) ^D,. (k) , (K 

J (27r) 

for the vertex function, and 

S (Pi) = Ql [ 4^,D^, (k) YS^ (pi - k) 7^ (59) 

J (27r) 

for the self-energy function. The fermion propagator is the usual form 

iS^ {x\ - x') = (0| T^{x[)^{x') |0) . (60) 

Its Fourier transform is defined by 

The two-photon exchange matrix element includes box and crossed-box matrix element 

-M':...... (p> q) = -^^°:..... (p> q) + MZ'Z.,, (P, q) , (62) 

where 

Ml:i^,^, (p, q) (63) 

— u (p. Si) (-iQir(^" (pi - A;)) u (q, ai) iD^^c^pV (-q, CTs) (^-iQ2rJ'' (P2 - k)j V (-p, Ss) , 

Ml;^Z^, (p, q) (64) 

— M (p. Si) (-iQir^" (pi - A;)) u (q, di) i/^/.^a/jV' (-q, as) (^-iQsr^ (^2 + k)j V (-p, Sa) . 

In the above we used the notation 

Df^uais = D^^, {k) {Dap {pi -qi- k) + Dp^ {P2 - q2 - k)) , (65) 

r^"(Pi-fc) = rS^ip,-k)^^, (66) 

r^;iP2-k) = 7%(p2-A;)7^ (67) 

rf{q2 + k) = YS^iq2 + k)^^. (68) 

The one-loop renormalization scheme in our formalism is considered in Appendix B. 

Note that the A^-matrix arises naturally in this formalism, i.e., Ai is not put in by hand, 
nor does its derivation require additional Fock-space terms in the variational trial state (41), 
as is the case in previous formulations (e.g. [14], [15]). 

10 



It is of interest to show that in the non-relativistic case the variational equation (47) 
reduces to the Schrodinger equation. Indeed, in the non-relativistic limit, the functions -Fsisj 
can be written as 



F,,,,(p) = F(p)A,, 



S2 5 



(69) 

where the non-zero elements of Aj^ for total spin singlet {S = 0) states are A12 = — A21 = 4j, 
while for the spin triplet {S = 1) states the non-zero elements are An = 1 for rris = +1, 



A 



A. 



-j= for nis = 0, and A22 = 1 for nis 



— 1. We use the notation that the 



12 = ^i.21 

subscripts 1 and 2 of A correspond to rris = 1/2 and m^ = —1/2 (or f and |) respectively. 
Substituting (69) into (47), multiplying the result by A^^^j and summing over si and S2, 
gives the eigenvalue equation, which determines the two-particle energy E 



{uj{p)+n{p)-E)F{p) 



{2nf 



rf3q/C(p,q)F(q), 



where 



/C(p,q) 



S1S20-10-2 



^s,sMsLa,a, (P, q) A.1C 



(70) 
(71) 



To lowest-order in (|p| , |q|) / (7721,7712) (i.e. in the non-relativistic limit), the kernel (71) 
reduces to /C = 51^2/ |p ~ ql ? and so (70) reduces to the (momentum-space) Schrodinger 
equation 



P^ 



2777^ 



e F(p) 



(27r)3 



rf^q- 



ri^(q) 



(72) 



|p-q| 

where e = E — {mi + 1712) and rrir = 77717772/ (^1 + ^2) is the reduced mass. This verifies 
that the relativistic two-fermion equation (47) has the required non-relativistic limit. 



5. Partial-wave decomposition and classification of states 

In the relativistic case we shall not complete the variational procedure in (47) at this 
stage to obtain final equations for the four functions -Fg^^j, because they are not independent 
in general. We require that the trial state must be an eigenstate of the relativistic total 
angular momentum operator, its projection, and parity, i.e.. 



J2 

J3 
V 



W trial) 



J{J+l) 

mj 
P 



W trial) 



(73) 



The total angular momentum operator is defined by the expression 



( rf^xT/;"^ {x) (L + 8)7/; (x) + j rf^x 0"f {x) (L + S)0 (x) 



(74) 



where L is the orbital angular momentum and S - the spin operator: L = x x p and 
S = I 0" . Using the field operators ip {x) and (x) in the form (18), (19), we obtain after 
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some calculations the expression for operator J. Explicit forms for the operators J^, J3 are 
given in papers [1], [2]. 

For a particle-antiparticle system charge conjugation invariance represents an additional 
requirement, i.e., we wish to construct the eigenvalues of 

C |e+e-) =C |e+e~). (75) 

However, this does not apply to the unequal mass case mi 7^ 1712- The functions -F^^sjIp) can 
be written in the general form 



where Y^ "^"^(p) are the usual spherical harmonics. Here and henceforth we will use the 



S2 



p 

notation p = |p| etc. (four- vectors will be written as p^). The orbital indices £^152^'^^ m^^ 
and the radial functions fs'iS2^ "^"'^ {p) depend on the spin variables Si and 82- Substitution 
of (76) into (41) and then into (73) leads to two categories of relations among the adjustable 
functions F^^sjIp)- 

Mixed-spin states 

In this case (!-siS2 = ^ = ^ and the general solution of the system (79) is 

F,,,,(p) = CiFS)(p) + C2FiS(p), (77) 

where Ci and C2 are arbitrary constants. Fslil{'p) and isis^p) are functions, which corre- 
spond to pure singlet states with total spin 5 = and triplet states with S = 1 respectively. 
The singlet functions have the form 

^Si(p) = C7(^^)---/^(p)y7»^»^(p). (78) 

where the Clebsch-Gordan (C-G) coefficients C^'9)'^^i-2 are: C(^9)™ii = C'^^a)^'^^ = 0, ^(^9)^12 = 
—(j{sg)m2i _ ]^_ Yox the quantum numbers mg^g^ one obtains: mn = 77222 = and mi2 = 
17121 = 'mj- The spin and radial variables separate in the sense that the factors fs^si (p) have 
a common radial function f"'{p). 

The triplet functions have the form 

fSAp) = C^S/™7'(P)>T""(P), (79) 

where Cj^ ^^ are the C-G coefficients for S* = 1, and 

"^11 =mj -1, mi2 = 7T721 = mj, 17122 = mj + 1. (80) 

We need to note that (78) is true for the singlet states J > 0, while (79) is true for the 
triplet states J > 1. Thus, the coefficient C2 in (77) is zero when J = 0. In other words, for 
J = 0, only the pure singlet state arises. For a system-like positronium the requirement (75) 
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decouples the singlet and triplet states for all J . Indeed, the charge conjugation eigenstates 



are 



\S9) = E^i™r"" jd'vf\p)Yr{msA^s, |0) 



with C = (— 1) for the pure singlet states, and 



\tr) 



E„{tr)ms^s 



d'pf{p)Yr{msA^sjo) 



(821 



with C = (—1) ^ for the pure triplet states, as it discussed in Appendix A. 

The states (81) and (82) diagonalize the Hamiltonian (13). Thus, for positronium-like 
systems, the states can be characterized by the spin quantum number S, and the mixed 
states (77) separate into singlet states (parastates S* = 0) and triplet states (orthostates 
S = 1). For distinct particles {mi ^ 7712) C is not conserved and there is no separation 
into pure singlet and triplet states in general. Thus for arbitrary mass ratio we need to 
diagonalize the expectation value of the QFT Hamiltonian AH = H — Hj^^ — M, in the 
basis of the states \sg), \tr) with (78) and (79) respectively for J 7^ 0. This can be achieved 
by the following linear transformation 



\s9q) 



U 



\sg) 

\tr) 



(83) 



where ?7 is a unimodular matrix Uu = U22 = a, f/12 = —U21 = — &, the components of which 
are, evidently, defined by the dynamics of the system. The new states, which diagonalize 
the expectation value of H, shall be called quasi-singlet | SQq) and quasi-triplet | tr-q) states 



\s9q) 



IN) 









d'vr\p)Yr{ms.Di^sM. 



d'pr\p)YT''{v)hlDl^sA^), 



M) 



where the coefficients C 



ySs jJ'lTis-^ S'2 



Jmj 



and C 






c 



Jnij 



{st)Jms^s2 
Jmj 

E 

v-iU2mj 



satisfy the following condition 

2 



C 



(st)Jmv 
Jmj 



2 (2 J 



^5) 



Note that these coefficients differ from C-G coefficients, due to the nature of the coupling. 
These coupled quasi-states arise only for J > 0. For J = pure S* = states occur. Quasi- 
singlet and quasi-triplet states are both characterized by the same quantum numbers J, mj 
and P = (—1)"'"'"^. Due to the unimodularity of the matrix U we can identify quasi-singlet 
and quasi-triplet states by quasi-spin (like isospin) t = 1/2 with t^ = =Fl/2, which is a 
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new quantum number. However, for our purpose it is more convenient to use the value 
s = ts + 1/2, which gives s = s^ = or s = St = 1 for quasi-singlet and quasi-triplet states 
respectively. In this case the labels Sg and St reflect better the meaning of the indicated 
quasi-states. As was shown in [2], for positronium the quasi-states become true singlet and 
triplet states with different charge conjugation quantum numbers. It is useful to note for 
subsequent calculations that the coefficients Ci and C2 in (77) are Ci = a, C2 = —b for 
quasi-singlet states s = Sg = 0, and Ci = b, C2 = a for quasi-triplet states s = st = 1. 

The triplet ^-mixing states 

These states occur for isiS2 = i = J ^ 1- The adjustable functions have the form 

Fs^sM = cf2'f-'^-^ f^\p)YT-r (p) + cfXf"-'^-^ f^\p)YT;r (p), (86) 

where the coefficients C/^ ^ are precisely the C-G coefficients. Expression (86) involves 
two radial functions f'^~^{p) and f'^'^^ip) which correspond to £ = J — 1 and £ = J+ 1. This 
reflects the fact that the orbital angular momentum is not conserved and l is not a good 
quantum number. The system in these states is characterized by J, mj^ and P = (— 1)"^. 
In spectroscopic notation '^^'^^Lj, these states are a mixture of ^ ( J — 1) j, and '^ ( J + 1) j 
states. The exception is the state with J = 0, for which the orbital angular momentum is 
conserved. Indeed, for J = the function f'^~^{p) does not exist (see Appendix A), thus 
the function Fs^s2{p) is defined only by the second term in (86). Note that £- mixing states 
appear for principal quantum number n > 3 only. 

For practical applications it is convenient to express pure states through the Dirac's F 
matrices, namely 

3 

Fs.s2 (P) = E /'(p) Vir- (P) V^-ps.. (87) 

The form of F^ (p) depends on the particular states and the index i corresponds to three 
cases: i = 1, when me = mj — 1, i = 2, when rrii = mj, and i = 3, when me = mj + 1. It is 
given below for the following cases: 

Pure Singlet States £ = J, J > 0, P = (-1)^+^ 

Ff (p) = F^ (p) = ^'Y-' (p) , (88) 

that is 

Wtr^al) = J^ /" ^'p /^(p)^?-' (p) Vi7 V_p,,6t^, D^p,^ |0) (89) 

SlS2 

Pure Triplet States (i = J, J > 0, P = (-l)^+i 
When me = mj — 1 

pjpx ^ \f mim2 \ ^'"^ A J + mj) ( J - mj + 1) ^ ^'"^ 

^^^^ ~ 2 V(cc'p + mi)(f]p + m2); V J{J+l) 

x(y+^7^-z(a°^+za°^))r7-Ap), (90) 
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when rrii = mj, 



'i rP) ^ - { ^^,,:^;i^,J "' ^TiTTTir ^-' ' '^"^ """ <'^' ' <"' 



and m^ = rrij + 1, 



17111712 



1/2 



(J-mj)(J + mj + l)V^' 



x(y-^f-.(a°i-.a°2))y7'^+^(p). (92) 

Pure Triplet States i = J - 1, J > 1, P = (-1)-^ 
When Tne = ttij — 1 

pJ_w-^ 1/^ mim2 \^'^ f {J + TTij - 1) {J + mj)^ ^'^ 



2 \(cjp + mi) (fip + ma)/ V J(2J-1) 

X (7^ + ^7^ - ^ (a°^ + ^a°^)) >T-r' (P) , (93) 

when rrif = mj 

' ^^^~ Vu, + m,){n, + m,)) \ J(2J-1) ) ^^ i^.-i(P)' 

(94) 
and TUi = TTij + 1 

pj-w-x _ 1/ 17111712 Y^"^ f {J - mj - 1) {J - mj)^ ^^^ 

^ 3 IPJ 



2 \{ujp + mi){np + m2)J V J(2J-1) 

x(7^-^7^-z(a°^-.a°^))y7_r^(p). (95) 

Pure Triplet States £= J+1, J > 0, P = (-1)^ 
When m£ = rrij — 1 

pJ+w^^ _ If mim2 y^^ /'(J-mj + l)(J-mj + 2)^ ^/^ 

■■^ 1 IPJ — 



2 V(tCp + mi)(f]p + m2)y V (J + l)(2J + 3) 

x(7^+.7^-^(a°i+^o-°^))F;--^(p), (96) 

when 772^ = mj 

pj+i .:^x / ^i"-^2 ^ ^'^ ( {J -mj + l){J + mj + l) \ ^'^ ,3 03. 

^2 (P) = 7 , wo I V / T I .WOT I QA (7 - «^ j >j+i (P) 

\(to'p + mi) (i2p + m2)/ V (J+l)(2J + 3) / ' 

(97) 
and 771^ = mj + 1 

^j+ip. _ 1/ mima V^' AJ + mj + 2) (J + mj + P ^ ^/' 

^^ ^"^^ 2 V(c^P + mi)(fip + m2)/ V (J+l)(2J + 3) 

X (7^ - .7^ - ^ (^°^ - ^^°^)) VT/i^^ (P) • (98) 
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6. The relativistic radial equations for two-fermion systems 

It is not possible to write one universal two-fermion radial wave equation, because the 
adjustable functions have different form for different states. Thus, it was important to 
classify all states of the system before deriving final radial equations. Now we return to the 
variational equation (47) from which we derive the radial equations for different states. We 
start with a particle-antiparticle system. 

6.1. Two-fermion wave equations for positronium-like systems 

It follows from the above analysis, that three sets of radial equations arise for this case. 
Singlet states £ = J, J >0, P = (—1) "*" , C = (—1) the radial equation is 



{2cu,-E)fip) 



m 



N{2 



TT 



q'^dq 



UJpUJq 



}Cip,q)f{q), 



(99) 



with the kernel 



/C (p, q) 



-'^^ Yl I dpd'^U^<T^^i (q) V-^a2Ms,s2a^a2 (P, Q) V-ps^^'/ (p) Ups„ 

i,j = l siS2cria2Tnj 



(100) 
and the normalization factor (we assume that the radial functions f'^{p) have been normal- 
ized) 



r ^ - 

N = I d^p 5^ 5^ wps.rf (p) \/_p,,\/_p,,r'/ (p) Up,, 

i=l siS2 



(101) 



3 

/ d'p J2 Tr 



,J /-N 7>a + "^-pJ (--\ ^"P" " ^ 



rr(p 



2m 



r/ (P) 



2m 



where F'^ is the matrix Fj of eq. (88), but with complex conjugate spherical functions, 
Y^i* (p). The four- vector p^ is defined as Pa = {Ep, — p). 
Triplet states i = J, J > 0, P = (-l)^+\ C = (-1)"^^^ 

For these states the radial equation formally coincides with (99), however the form of the 
F/ matrix must be taken from (90)-(92). 
Triplet states i= J±l, P= (-1)*^, C = (-1)^ 

In this case the variational equation (47) leads to a system of coupled equations for the 
two independent radial functions f"'~^{p) and f'^~^^{p)'- 



{2cOp -E)¥ (p) 



m 



(27r) 



q^dq 



UJpUJq 



K(p,g)F(g), 



(102) 



where 



F(p) 



N-^-^f-^-\p) 
NJ+^fJ+\p) 



(103) 
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and 



The kernels /Cm„ are 



K{p,q) 



/Cii (p, q) /Ci2 (j5, q) 

^21 (p, g) /C22 (p, g) 



(104) 



^mn (P,g) = -^ X] Yl / C^PC^q^q^iri" (q) t;-qa2 A<siS2<Tia2 (P, q) ^-psar^^" {p) Ups^, 

i,j=l (Ji(J2SiS2mj 

(105) 
Here ii = J — 1, £2 = ^ + 1- The T^" matrices are defined by (93)- (98). The system (102) 
reduces to a single equation for J = since f'^^^{p) = in that case. The normalization 
constants A^"'^^ are defined by (101) with corresponding matrices T- ^ (99)-(104). 

6.2. Two-fermion wave equations for muonium-like systems 

After completing the variational procedure we obtain the following results: 
For the states with £=J = 0, P = — 1 the radial equation is 

17111712 f q^dq 



iujp + n,-E)f'ip) 



N{2tt) 



CUpCUqi lp\ Iq 



ICip,q)fiq), 



(106) 



where the kernel /C {p, q) is defined by 

^ (P, g) = -^ Y j dpdqu^^^-f^ (q) v_^^^Ms,s2aja2 (P> q) 



v-ps2l^ip)upsi- (107) 



S1S20-10-2 

The normalization factor is 



N = — I d^pTr 



5 7>al + ^l 5 7>a2 " "^2 

7 ^ 7 



:io8) 



2mi 2m2 

For quasi-singlet and quasi-triplet states (J > 1) we have the system of two equations 

q^dq 



where 



and 



ujp + np-E)¥ (p) = 

¥ip) 

K(p,g) = 



17117)72 



(27r) J y^UpUqflpQ, 



]\[isg)fisg)Jtp\ 

Ar(*0/(*'-)-^(p) 



}Cn{p,q) ICi2{p,q) 
/C2i(p,g) /C22(p,g) 



:K(p,g)F(g), 



(109) 
(110) 
(111) 



The kernels JCmn (p, q) are 



^mn {P, q) = -i J^ / d^pd^q Mq^iT" (q) V_^„2Ms^S2a^a2 (P, q) ^-psar"" (p) Mpsj • 

cricr2SiS2mj 

(112) 
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Here we used the following notation 

T'ip) = r^(p), r'l (p) = r'^ (p) , 

3 3 



;ii3) 



i=l 



i=l 



The normalization factors are 



j^isgUtr)^ d^^Y^Tr 



i=l 



^,j ^^^ rP.i+m,^j ^j5) rp.2 - m. 



2mi 



2mi 



(114) 



For the triplet states with i = J^l, we have two independent radial functions /"^ ^{p) and f'^^^{p)- 
Thus the variational equation (53) leads to a system of coupled equations for /'^"^(p) and f'^^^ip) 



where 



and 



(a;p + fip - E) F (p) = 
¥{p) 
K(p,g) = 



17111712 



(2vr) 



q^dq 



UJpUJ q\ lp\ Iq 



K(p,g)F(g) 



NJ-^f-^-\p) 



/Cii (p, g) /Ci2 (p, g) 
^21 (p, g) /C22 (p, g) 



(115) 
(116) 
(117) 



The kernels /C^n are 



/c„„ (p, g) = -^ Xl m / #c?qMq<,,r." (q) V-^„^Ms^s2<T^a2 (p, q) ^-psar^^™ (p) ^psi, 

i,j=l cr\cr2SiS2mj 

(118) 
where £i = J— 1, £2 = ^+1- The normalization factors are defined by analogy with (114). 
The system f)115p reduces to a single equation for f^^^ijp) when J = 0, since f'^~^{p) = in 
that case. 
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7. Radiative corrections to O (a^) for arbitrary mass ratio 

The relativistic radial equations derived above cannot be solved analytically, so an ap- 
proximation method (numerical, variational, or perturbative) must be used. We shall use 
first-order perturbation theory to calculate O (a^) corrections to the energy for some states. 

As mentioned above, at the one-loop level radiative corrections include vacuum polar- 
ization, fermion self-energy and vertex corrections. We calculate their contribution to the 
energy shift of singlet states with i = J = 0. These states are described by the radial 
equations (99) and (106) for positronium- and muonium-like systems respectively. Since the 
solution of (106) reduces to that of (99) for two equal masses, we concentrate on equation 
(106) only. The energy eigenvalues En^j can be calculated from the equation 

/»oo /»oo 

En,j dpp'f{p)f\p) = / dpp'{oo, + nj,)f\p)f\p) 
Jo Jo 

m.m, r PV'^f'^ .X(p.,)f(p)f(,). (119) 



N (27r) Jo ^JuJpVLpUJqVLq 

which follows from (106). For radiative corrections to O (a^) we take f'^{p) to be the non- 
relativistic hydrogen wave functions in (119) and obtain the result 

2 

Ae(a^) = E„^, - (^1 + ^2) + ^ - A5 (a^) (120) 

m,m, ^ r PVdpdq (/c--(p,g) + /C--(p,g) + /C^ -(p,g)) /^(p)/'^(g). 



N{2'!t) .tl^2>^0 ^/idp^Ldq^q 

where rrir = 17111712/ {rrii + 1122) is the reduced mass. As (a^) are energy corrections to O (a^) 
derived in [2]. The kernels in (120) correspond to the radially reduced form (107) of the 
second order matrix elements (53)-(56). They are, explicitly (/i^ = a/Airmi): 

^_^ QiQ^ /-^^^^^ (121) 

647rm^m2 J 
xTr [{-f^pix + mi) A'^ {pi - q^) (7^^^ + mi) 7^ {-f%x - "^2) Y {l%x - "^2) 7^ , 



^'""''{P.q) = Z'^t\ I dvdv,D^,{p-q)W'^{p,-q,) (122) 

b47rm^m2 J 

xTr [{-f^pix + mi) Y {l^Qix + "^1) 7^ (7^g2A - "22) 7" (7^P2A - "^2) 7^ , 

1 2 •-' 

xTr [(7^PiA + "^1) S (pi) S^ (pi) Y (7^giA + mi) 7^ (7^g2A - ^2) Y {l%x - "^2) 7^] • 

The contribution of the kernels JC^'^^^, /C^"^2 jg obvious, and follows from the symmetry 
with respect to mi and m2. As will be shown, the contribution of /C™"**i, /(;™assi ^ j^mass^^ 
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and /C^^^^a is zero. The renormalization of the functions A'' {pi — qi) (64), 11'^" (pi — gi) 
(57), and S (pi) (59) imphes that they should be replaced by the following expressions (the 
justification is discussed in Ref. [16]), 

^' iPi - qi) ^ ^^^^'^ (Pi - ?i). ' (124) 

n^" (Pi - gi) - -^^"y^^ (Pi - gi)\ (125) 

S (pi) ^ 0. (126) 

After that the calculations are straightforward and we obtain the non-vanishing kernel con- 
tributions 

IC ' (P, q) = -a — ^—^ ^5j,o, 127 

^^""(P,g) = ^^^,o. (128) 

To calculate energy corrections we use the nonrelativistic hydrogenic momentum-space radial 
wave function f'^ip) (ref. [17], eq.l25). The corrections are given by 

a^mr 2 (2mi - m^) m^ 

^^ = — § ; ^2—^^,0, (129) 

n-^ 71 (mi + 1712) 



^^'"'^^ = -^TTZ^^^J,o, (130) 



Aa^rrir ml 
15?™'^ ml 

These results agree with Ref. [18]. 

It should be noted that a similar treatment of positronium and muonium has been con- 
sidered by Zhang and Koniuk [16, 19]. They used postulated equations with an inserted 
invariant Ai matrix. These authors show that the inclusion of single-loop diagrams yields 
positronium energy eigenstates which are accurate to 0{a^,a^ In a). In the present varia- 
tional treatment the equations and results are derived from first principles. 

8. Concluding remarks 

We have shown that the variational method can be formulated in a way that allows one 
to derive relativistic few-body equations, which can include interactions to any order of the 
coupling constant. The method is based on a reformulation of QED, in which covariant 
Green's functions are used to solve partially the underlying Euler-Lagrange equations of mo- 
tion. This leads to a Hamiltonian which contains the Green's function sandwiched between 
fermion currents directly (equations (14)-(16)). The eigenvalue equation P^ {ip) = P^ {ip), 

where P'^ = iH,P] is the energy momentum operator of the QFT is formulated variation- 
ally, 5 ({tpl H — E |'?/')j=jj| ] = 0. Time evolution operators are used to recast the problem in 
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such a way that Feynman diagrams of any order can arise in the kernels that describe the 
interparticle interaction (eq.(39)). 

We illustrate the utility of this formulation by deriving relativistic momentum-space 
wave equations for two-fermion systems like muonium and positronium. These equations 
(47), which become (99), (102) for particle-antiparticle systems, and (106), (109), (115) for 
muonium like systems, describe the behavior of the two-fermion systems in principle to all 
orders of the coupling constant for arbitrary mass ratio. For bound states of the two-fermion 
system the trial states are chosen to be eigenstates of the total angular momentum operators 
J^, J3 and parity, and also of charge conjugation for particle-antiparticle systems. A general 
relativistic reduction of the wave equations to radial form is given for arbitrary masses of the 
two fermions. For given J there is a single radial equation for total spin zero singlet states 
(106), but for other states there are, in general, coupled equations, (109) for mixed-spin 
states, and (115) for triplet ^-mixing states. We have shown how the classification of the 
states follows naturally from the system of eigenvalue equations (73), given our trial state 
(41). 

We use the derived radial equations to obtain approximate perturbative solutions for 
the two-fermion binding energy to O (a^) in the fine-structure constant for all singlet states 
with total angular momentum quantum number J = £ = 0. Results for other states can be 
obtained in an analogous manner. 

The method presented here can be generalized for systems of three or more fermions. 
This shall be the subject of a forthcoming work. 
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Appendix A: Matrix elements for high-order corrections 

To classify the corrections we write out the Hamiltonian density Tin {x) in explicit form 



rtii ix) = Hr, (x) + Hi^ (x) + Hr,,^ (x) + Hi^,, (x 



(p^ y^ / 1 



where 



Ql 



'^li, (^) = -y" / d'^x ^{x)-i^^{x')D^^{x - x')tlj{x)Yi^ix) 



Ql 



T-Cr,,^ (x 



i"f> 



Hr^,:. (x) 



^tpip 



^^^^ ' d^'x' ^{x'h>'^{x')D^,{x - x')(f){x)Y<l){x) 



' d!^x 4){x)'^^(t){x)D^j,{x — x')ip{x)'y^'ip{x) 



(131) 

(132) 
(133) 
(134) 
(135) 



First, we consider the element T ( T-Ci (xi) Ti-i^, (x) 5 (t)). It is not difficult to show that 



i^triail / d^xid^xT [Hi. {xi) Hi., {x) 6 (t) ) l^pt^^J 



Q\Q2 



(136) 



d xid x'ld x'd xD^i, {xi — x'l) Daf3{x — x)5 (t) 



x2(^ 



trial I 



/ : ip{x'i)'^^ip{x'i)Tr {iS.^ {x' — Xi) •y'^iS^ {xi — x') 7") (f){x)'y^(f){x) : \ 
+ : ip^XijYiS^ {xi — x') '~i"'iS^ ix' — x\) 'y'^ip{x[)(f){x)'y^(f){x 
■i/j{x[)'y^iS^ {x'l — Xi) •y^iS^ {xi — x') 7"-?/'(x')0(a;)7^0(x 
ilj{x')y^iS^ {x' — xi) y^iS^ {xi — x\) y^ilj{x'i)(f){x)y^ (f){x 



+ 

V + 



W trial) 1 



where we used the symmetry of the Green's function Df^^ (xi — x'^) = Di,f^ {x[ — xi) and the 
standard expression for the fermion propagator (60). Using the Fourier transform (18), (19), 
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(49), (61) we obtain 



(^iria/l / d^Xid'^xT (n^ (Xi) ili^^ {x) 5 (t) j l^trial) 



(137) 



QIQ2 mim2 



(Ppi(fp2d^Pid^p'2 f d'^k 



2 (2r^6 ^ 



X 



in 



X 



X 



X 



X^sls^ (Pl' P2) ^siS2 (Pl> P2) ^ (Pl - Pi + P2 - P2' 

/ Df,^ (p; - Pl) Dafs (K - Pl) X _ \ 

xu (p'l, s',) 7^1. (pi, si) Tr [5^ (A; - p[) Y S^ {k - p^) 7"] V (p2, S2) l^V (p'^, 4) x 

+D^^ (A;) Dap {p'l - Pl) X _ 
u (p'l, s;) 7^^^^ (A; + p[) rS^ {k + Pl) Yu (pi, siV^^ (p2, S2) i^V (p'^, 4) 

+11 (p;, ^i) rs^ {p'l - k) rs^ (p'l) ru (pi, si)v_ (p2, .2) i^v (p'^, 4) , 

V V +m(p;,s'i)7"^^(Pi)7"^v(^ + Pi)7^w(Pi,Si)\/(P2,S2)7''V^(P2,4) / / 

The calculation of matrix elements is presented here in an arbitrary frame. In the rest frame 
the adjustable functions Fs^s2 (Pi; P2) reduce to Fs^s2 (Pi) ^ (Pi + P2)- 

The next term Tii^ {xi)'Hi^^ (x) yields a similar result, but with the Green's function 
Dai3 {p'l - Pl) replaced by D^p {p'^ - ^2)- Finally, the terms Hi^ (xi) Hi^^ {x)+Hi^ (xi) Hi 
and Hi^^ (xi) Hi^ (x) + Hi^^ (xi) Hi^ (x) give the same contribution and can be combined 
together in the form 



X 



Vrtrial 
17111712 






Z 



(2>r)' ,^ 

^ ' S\s'2SiS2 

A2)vaci 



d^Pid^P2d^p'id^p'2 , , ,. , . , , i\t ■\ 
^^44 (Pi' P2) ^^1^2 (Pi, P2) 5 (Pi - Pl + P2 - P2) (-«) X 

(^Pl ^p'^ ^P2 ^V'2 ) 



,(2)veri 



-^is;k (pi> p'i> P2, P2) + -Mir-, (Pl, p'l, P2, P'2) 

, , ,(2)mass'-. , , i \ . a A2)mass'J , , ,\ 

+-^s[s',s,S2 (Pl' Pl' P2, P2) + ^s[s',s,S2 (Pl' Pl' P2, P2) 



138) 



j(2)vaci 



, {2)veri 



, (2)masSj . ,{2)mass" 



where M'rr' , A^^T^r^ , A<77""S MVr'^' are defined by (53)-(56). 

S\S2S1S2^ s\S2SlS2^ s\S2SiS2 ' s'j^S2SlS2 •' ^ ^ ^ ' 

Similarly we can show that the next terms contribute to the two-photon exchange process 



vPtrial 

qiq2Tnim2 
2 (27r)^ 



d'^xid'^x T ( ( '^i'*^ ^^^'^ '^i'^^ ^^'^ ^ '^i'^^ ^^^'^ '^i^'^ ^^'^ 
+^/^^ (xi) iii^^ (x) + H/, , (xi) Hi., (x) 



E 



(i^Pi(i^P2(i^Pi(i^P2 



t^2 r / I / / \ 

XT^o (p2 - Pl + Pl - P2) 



S'^S-2S1S2' (^p'i^Pi^P2^pO 
X^sls'2 (Pl, P2) ^s^S2 (Pl, P2) i-i) -M^s^sisa (Pl, Pl, P2, P2) 

where .Mj^'g g. is defined by (63)- (64). 



Ht)] \^tr^al) 
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Appendix B: The one-loop renormalization scheme 

It is well known that the vacuum polarization function (57) can be written as 

n^" {p[ - pi) = [{p[ - p,r ip'i - PiT - ip'i - Pif g"""] n ip[ - p,) , 

where 



U{p[-p,) = U{0) + U{p[-p,) 



ni(o) 
n(p'i-pi) 



127r4 
27r2 I 



d^k 



{k? — iTii + ie 
dzz (1 — 2;) In 



>2' 



1- 



{p'l-Pif z{l- z) 



2 

mi 



The infinite sum 



iD'^p{p'l-Pl^P'2-P2 



(139) 

(140) 
(141) 
(142) 

(143) 



D^t3 {p'l - Pl,p'2 - P2) + D^y {p\ - Pi) W {p\ - Pi) iDat3 {p'l - Pl,p'2 - P2] 

+D,, {p[ - Pi) W" {p[ - Pi) D^s {p'l - Pi) n'^ {p'l - Pi) iD^p {p'l - pi,p'2 - P2) 

gives the dressed Green's function, if the vacuum polarization is due to the first particle. 
The bare Green's function D^^ (p'^^ — pi,P2 — P2) is defined by (50) and, in Feynman gauge, 
has the form 



D^p {p'l -Pi,p'2-P2) = g^p 



(144) 



{p'l-pif -{p'2-P2y 

Dropping terms like {p'^ ^ Pi)"^ {p'l — Pi)'^ (they contribute zero when they are contracted 
into the final or initial currents) we obtain 

iD'^p{p'i-Pi.p'2'P2) (145) 

^ (1 - n {p\ - pi) + n^ {p\ - pi) - n^ {p\ - pi) + ..) iD^p {p[ - pi,p^ - P2) 

iDf,fs{p'i -Pi,p'2-P2) 



l + U{p[-pi] 
If we take into account polarization due to the second particle we obtain 

iD'^(3{p'i-Pi^P'2-P2] 



(146) 



l + U{p' 



1 1 

— -iD^p {p[ - pi,p'2 - P2) ^ , „. , 

.Pi -pi) 1 + n {p2 



'2 - P2) 



= iDnR (p\ — Pi, P'r, — P2) = 

i + ni(o) + n(K-Pi) '^^' ' "^'^ i + Yi2{^) + n{p',-p2] 
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Then we define 

^ l + ni(0), -1- = 1 + I[2i0), (147) 



7(1) ' ^^^V"^' 7(2) 



and obtain 



^-C'ufl (Pl - Pi , P2 - P2) = ^mJ:, , :^^M/3 (Pl - Pi > P2 - P2^ ^ 



(148) 
Introducing the renormalized quantities 

ni?(p'i-pi) = zi'^Ti{p[-p^), (149) 

n^ (P2 - P2) = 4^^n (p2 - P2) , 



expression (160) can be written as 

j+^R {p[ - p,)) (1 + Ur ip', - P2)) 



7(1)7(2) 

^3 ^3 



^^)./3 (P'l - Pi, P2 - P2) = r, , TT /^/ ^ \M. , TT /^/ Z:\\^^I^P (Pl - Pi' ^2 - P2) 



(150) 
mahzation of charges. In the lowest 
of expansion, the dressed Green's function (150) is 

)' 
Note, that in the nonrelativistic limit the IIr- functions take on the form 



Constants Z^ and Z3 shall be absorbed by renormalization of charges. In the lowest order 



iD'a (Pi - Pi> P2 - P2) = Z'i^zf^ (1 - n^? {p\ - pi) - Ur {p2 - P2)) iD^p {p[ - pi.p'^ - P2) 



(151) 



TT / / X «(Pl-Pl)^ TT ^ ' ^ "(P2-P2)^ i,^^. 

Next, we consider the vertex function 

r"(p'i,Pi)=7" + A"(p'i,Pi), (153) 

where A" can be written in the form 

A" (p;,pi) = ^1 (p;,pi) 7" + ^2 (p;,Pi) """^ ^^^~^^^^ (154) 

where the scalar function J-'i (pi,pi) diverges, while JF2 (pi,pi) is finite. Therefor, 

r"(p'i,pi) (155) 



io-"^ (p\ — Pl)p, 

7" + ^1 (p ,pi) 7" + ^2 (p ,pi) — r ^ 

(i + ^i(p;,p;))7" + (^i(p'i,Pi)-^i(p'i,p'i))7" + -^2(p'i,Pi 



ia"^(p;-pi)^ 



2mi 
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Then we define a value 

^ = i + ^i(p;,p;), (156) 

which will be absorbed by the renormalized charge 

Qm = Qi^ (157) 

The renormalization of the self-energy can be treated similarly. 
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